Tuning the quantum oscillations of surface Dirac electrons in the topological insulator 

Bi 2 Te 2 Se by liquid gating. 
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In Bi2Te2Se, the period of quantum oscillations arising from Dirac surface electrons can be in- 
creased 6-fold using ionic liquid gating. At large gate voltages, the Fermi energy reaches the N = 
1 Landau level in a 14-Tesla field. This enables the i-shift predicted for the Dirac spectrum to be 
measured accurately. A surprising result is that liquid gating strongly enhances the surface mobility. 
By analyzing the Hall conductivity, we show that the enhancement occurs on only one surface. 
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A Topological Insulator (TI) is characterized by the 
existence of current-carrying surface states that traverse 
the bulk energy gap [l|, [|j • There is strong interest in the 
helical nature of the surface states, which results from 
the locking of the electron's spin transverse to its mo- 
mentum. In the bismuth-based TI materials, photoemis- 
sion spectroscopy [3[ and scanning tunneling microscopy 
(STM) [i[ have confirmed the spin-locking feature. In 
transport experiments, the surface states have been de- 
tected by surface Shubnikov de Haas (SdH) oscillations in 
Bi 2 Te3 [1] and (Bi,Sb)2Se 3 Q, and by Aharonov-Bohm 
oscillations in Bi2Se3 nanowires Q- 

Among the Bi-based TI materials, Bi2Te2Se currently 
displays the highest bulk resistivities (p = 1 — 6 Slcm at 
4 K) [8Hll|. Despite the large p, SdH oscillations may be 
tracked to temperatures T as high as 38 K @ . A persis- 
tent problem, however, is that the surface Fermi energy 
Ep in as-grown crystals is still quite high (~200 meV 
above the Dirac Point). An in-situ method that demon- 
strably tunes Ep would greatly facilitate experiments at 
the Dirac Point, as proposed in Refs. [11,11]. Several 
groups have applied conventional electrostatic gating to 
tune the chemical potential p in exfoliated crystals |12| - 
ITij and in thin-film samples of Bi2Se3. The newer 
technique of liquid gating has also been used on Bi-based 
materials [Ra-|20l]. However, in these experiments, SdH 
oscillations were either not detected at all or poorly re- 
solved. Tuning of the quantum oscillations and showing 
that they arise from surface Dirac electrons remain to be 
established. 

Here we report that the surface SdH oscillations in 
Bi2Te2Se can be tuned over a broad range using the ionic 
liquid (DEME-TFSI). In lowering E F substantially, we 
access the N = 1 Landau level in a magnetic field B — 14 
T. This allows the i-shift characteristic of Dirac electrons 
to be measured with greatly improved resolution. We 
find that liquid gating leads to strong enhancement of the 
mobility p s of the surface carriers. Aside from moving 
Ep closer to the Dirac Point, the tunability yields direct 
information on the surface and bulk conduction. 

As reported earlier 0, [t| [ll|, the resistance R(T) in 
Bi 2 Te2Se rises monotonically to very large values as T —} 
4 K (curve at Vg = in Fig. QJ,). Analysis of the Hall 



coefficient Rr at 5 K (Fig. [Tb) reveals a population of 
bulk n-type carriers much higher than the population of 
surface electrons. Nonetheless, a modest, negative gate 
voltage Vg can increase R by 40% (Fig. [TJi) and \Rh\ 
by a factor of 2 at 5 K (Panel b). Vq is applied to the 
gold electrode (inset in Fig. QJi) at 220 K, and the sample 
is then cooled below the liquid's glass transition. After 
the low-T measurements are completed, the sample is 
warmed to 220 K (at 2 K/min.) and Vg is reset. [2l[ 
At 4 K, the large E'-field induced by the surface anion 
density Ni on (l-4xl0 14 cm~ 2 ) creates a depletion layer 
that penetrates deep into the bulk (5-20 /mi). [23[ As 
shown in Fig. [T|b (inset), the induced upward bending of 
the bands decreases Ep. 

At each value of Vg, the curves of R vs. B display 
SdH oscillations. To focus on the SdH signal, we have 
subtracted off a smooth background ps to isolate the os- 
cillatory part of the resistance, Ap xx = p xx ~ ps- Figure 
[2^, displays plots of Ap xx in Sample 1 versus 1/B for 
5 values of Vg- The period of the SdH oscillations in- 
creases monotonically as Vg changes from to -4.2 V, in 
accordance with our expectation that Ep is decreasing. 
Surprisingly, the SdH amplitude is strongly enhanced be- 
tween Vg = and -2.1 V (the former is shown ampli- 
fied by 5x). The dotted curves are the best fits [llj 
to the Lifshitz-Kosevich (LK) expression for SdH oscil- 
lations using only one frequency component. From the 
fits, we may infer how the surface mobility p s changes 
with Vg (see below). The same trends are evident in 
Sample 2, which has a higher starting surface density n s 
but is taken to B = 45 T (Fig. [2p). We find that the 
SdH oscillations are not resolved at Vg — 0, but become 
prominent at Vq = -1-5 V. 

In finite B, the surface electronic states are quan- 
tized into Landau levels (LLs) with quantum numbers 
N = 0, 1, • • ■ . The index field B n is the field at which 
Ep falls between two LLs. For Schrodinger states, the 
integer n counts the number of occupied LLs (the highest 
filled LL has N max = n — 1). Using the level degeneracy 
Be/h per spin, we then have 1/B n = ne/(hn s ) (n s is the 
surface density, e the elemental charge and h is Planck's 
constant). 

For Dirac electrons, however, we have n + h filled LLs 
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FIG. 1: (color online) The resistance R per square (Panel a) 
and Hall coefficient Rh vs. T (Panel b) in Bi2Te2Se at se- 
lected Vg in Sample 1. Rh is measured at fixed B (3 T). 
Changing Vg from to -2.8 V increases R by 40% and \Rh\ 
by 2 x . The inset (Panel a) shows the cell housing the sample 
and the ionic liquid DEME-TFSI. The Au electrode (a circu- 
lar plate of radius 1.5 mm) is separated by 0.5 mm from the 
sample. [2l|] The inset in (b) is a sketch of the band bending 
induced by liquid gating. Negative ions deposited on the crys- 
tal leads to upward band-bending. At the surface, this causes 
Ef to decrease towards the Dirac Point. LLs are shown as 
solid half-ovals. 



when B = B n (now N max — n). The additional i derives 
from the N — LL, or equivalently, from the 7r-Berry 
phase intrinsic to each Dirac cone [24| . The relation be- 
tween 1/B n and n is now 1/B n = (n + ^)(e/hn s ) - a 
straight line that intercepts the n-axis at n = — h. In 
both cases, G xx is a local minimum at B n . 

If resistivity curves are used, B n should be identified 
with the maxima in Ap xx . This point is discussed in 
Refs. [ll|, [25] In Fig. [3k., we plot as solid symbols B n 
in Sample 1 against the integers n (the open symbols 
corresponding to the minima are plotted against n + ^ ) . 

At each Vg, the slope of the straight lines yields the FS 
area Sf- As \Vg\ increases from to 2.8 V, the slopes of 
the best-fit lines decrease by a factor of 6.4, reflecting a 
steep decrease in Sf- This decrease saturates when \Vq\ 
exceeds 2.1 V. 

In Panel (b) , we show the high- field behavior for | Vq \ > 
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FIG. 2: (color online) Traces of SdH oscillations in the re- 
sistance versus 1/B, showing systematic changes to the oscil- 
lation amplitude and period with gate voltage (bold curves, 
displaced vertically for clarity) . The dashed curves are fits to 
the LK expression with one frequency component [ill]. Panel 
(a) shows traces of Ap xx vs. 1/B at 5 K measured to 14 T for 
5 values of Vg (Sample 1) [22l ]. The largest increase in ampli- 
tude occurs between Vg = to -1.4 V. The curve at Vg — is 
shown amplified 5 x . All other curves share the same vertical 
scale. Panel (b) displays traces of AR XX vs. 1/B at 0.3 K 
measured to 45 T at Vg as indicated (Sample 2) [22| . Arrows 
indicate n = i (Ef at center of broadened N — 1 LL). 



2.1 V in expanded scale. At these large bias values, the 
intercepts cluster around n = — i (-0.46, -0.56 and -0.61). 
In the corresponding plots for Sample 2 (Panel c), we 
find that Sf decreases by a factor of 2 between Vg — 
-1.4 and -6 V. In the limit 1/B —> 0, the intercepts are 
at n = —0.35, -0.40 and -0.42. In both samples, the 
last feature observed at the highest B (minima in R xx ) 
corresponds to n m %n — \ (this implies that Ep lies in the 
middle of the broadened N — 1 Dirac LL). With such 
a small n m i n , we may rigorously exclude an intercept 
at n = in the limit 1/B ->• 0. HU Thus the index 
plots provide rather conclusive evidence that the SdH 
oscillations arise from surface Dirac electrons. 
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FIG. 3: (color online) Index plots of the integer n vs. 1/B n 
at selected Va in Sample 1 (Panels a and b) and 2 (Panel c). 
Maxima of Ap xx (solid symbols), corresponding to the index 
fields B n , are plotted against n. Minima (open symbols) are 
plotted against n + |. Panel (a): As Va changes from to 
-2.1 V, the slope of the best fit lines decreases 6-fold. Further 
increase in \Va\ leads to saturation. In Panel (b), the high- 
bias curves are displayed in expanded vertical scale. In the 
limit 1/B — > 0, the best-fit lines have intercepts at -0.46, -0.56 
and -0.61, consistent with Dirac electrons. The intercepts 
for Sample 2 (Panel c) also cluster near -0.45 in the limit 
1/B ->■ 0. 
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FIG. 4: (color online) Panel (a): The observed Hall conduc- 
tivity a xy vs. B in Sample 1, showing weak-B curvature at 3 
values of Va (curves displaced for clarity). At each Va, the 
outer curves are the data (solid black curve) and the fit to 
Eq. Q] (superposed blue dashed curve). The inner (red, solid) 
curve is the surface term G s xy /t fixed by n s and /j, s - (Weak 
SdH oscillations are apparent in the observed a xy .) The dif- 
ference between the outer and inner curves is the bulk term 
a xy . At Va = -4.2 V, G xy /t accounts for 83% of a xy in weak 
B. Panel (b) shows that, with increased gating, ^t s increases 
from 720 to 2,480 cm 2 /Vs while fi b stays very small (20-30 
cm 2 /Vs). Panel (c) compares the sharp decrease in n s with 
the mild change in rib with gating, 
saturates. 



When \V G \ > 2 V, 



Returning to Fig. [3J we have fitted the SdH oscillations 
to the LK expressions (shown as dashed curves). The 
damping of the oscillations versus B yields the surface 
mobility /i s . As shown in Fig. |4Jd, /i s in Sample 1 rises 
from 720 to 2,480 cm 2 /Vs as \V G \ is increased to 4.2 
V. [22| In Fig. HJ;, the decrease and eventual saturation 
in Sf is plotted as a surface density n s = k F /(4ir) (per 
spin). The saturation at large \V~g\ either arises from 
induced chemical reaction or from Ep meeting the top of 
the valence band. 26] 

One figure-of-merit in TI crystals is the ratio of the 
surface to bulk conductances r) = G s /G b in zero B (with 
G r = G^, x (0), r = s,b). In Sample 1, 77 ~ 0.05 is quite 



small (compared with 77 ~ 1 obtained in Ref. [9]). How- 
ever, in the Hall channel, the ratio t}h = G s xy /G h xy of the 
surface and bulk conductances (G xy and G xy , respec- 
tively) is enhanced by /i s //jfc, which can be very large. 
We define rib and fit, to be the bulk electron density and 
mobility, respectively, averaged over the whole crystal. 

As shown in Fig. a distinctive feature of a xy at low 
T is the curvature in weak B, which grows with increasing 
\V~g\- We may use the semiclassical 2-band expression for 

a xy = n s en s + + n b ert b B, (1) 
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where the first term is G xy /t, with t the thickness (50 /jm 
in Sample 1). With n s and fi s fixed by analysis of the 
SdH oscillations, this term is non-adjustable. The second 
term is the bulk Hall conductivity a xy in the low-mobility 
limit. With the sole adjustable parameter Pb = n^/i 2 , we 
find that Eq. [T] gives a very good fit (dashed curves). 
For comparison, we have also plotted G*/t (inner, faint 
solid curves). Combining Pb with the zero-B value of a xx , 
we finally obtain rib and fib separately for each value of 
Vg- These are reported in Figs. HJd and Hfc. The small 
values of fib (20-30 cm 2 /Vs) result in a large fi s /fib ~ 
100 and rfn ~ 5. This accounts for the pronounced low- 
B curvatures seen in Fig. |4] 

The analysis implies high-mobility Dirac electrons in 
parallel with a much larger population of bulk electrons. 
Because of the 100-fold difference in mobilities, the Dirac 
electrons produce 83% of the total weak-f? Hall conduc- 
tance at large \Vq\. The fits include the surface Hall con- 
ductance from only one surface. Since its G xy already 
accounts for most of the observed a xy , there is very little 
room left for a second surface term. We estimate that the 
Hall contribution from the other surface is less than 2% 
of a xy , which implies that its fi s is <300 cm 2 /Vs. This 
cannot produce resolvable SdH oscillations. 

The large enhancement of fi s by liquid gating (Fig. 0b) 



is perhaps the most intriguing feature reported here. To 
our knowledge, this is the first realization of enhancement 
of surface SdH amplitudes by an in situ technique. [13] A 
recent STM experiment [H| reveals that the Dirac Point 
closely follows spatial fluctuations of the local potential 
on length scales of 30-50 run. This could lead to strong 
scattering of surface electrons. We speculate that, under 
liquid gating, the anions accumulate at local maxima in 
the potential, thereby levelling out the strongest spatial 
fluctuations. [29| The results provide encouragement that 
alternative routes that even out local potential fluctua- 
tions can lead to further improvements in fi s . 
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